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Abstract The structure of a society depends, to some extent, on the incentives of the individuals they 
are composed of. We study a stylized model of this interplay, that suggests that the more individuals 
aim at climbing the social hierarchy, the more societys hierarchy gets strong. Such a dependence is 
sharp, in the sense that a persistent hierarchical order emerges abruptly when the preference for social 
status gets larger than a threshold. This phase transition has its origin in the fact that the presence of 
a well defined hierarchy allow agents to climb it, thus reinforcing it, whereas in a "disordered" society 
it is harder for agents to find out whom they should connect to in order to become more central. 
Interestingly, a social order emerges when agents strive harder to climb society and it results in a state 
of reduced social mobility, as a consequence of ergodicity breaking, where climbing is more difficult. 



1 Introduction 

The emergence of social elites has interested social scientists ever since Pareto's observation of persistent 
inequalities in our societies |14| . Inequality is acceptable if it results from differences of individuals in 
terms of their capabilities, but not if it results, in one way or another, from discriminatiorj^ Not only 
discrimination conflicts with ethical principles that all individuals are a priori equal and should have 
access to the same opportunities. It also damages societies in terms of efficiency (18j as it hampers 
social mobility, preventing society to promote individuals to positions in the social hierarchy, that are 
consistent with their capabilities. 

We introduce the social climbing game, a highly stylized model of a society, where individuals 
attempt to optimize their position in the network, by becoming as central as possible. The assumptions 
of the model are rooted on empirical and theoretical evidence coming from the social sciences. There, in 
the early years of network analysis, it was found that the importance of an individual within a network 
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is related to some quantification of liow central pil22) tfiis agent is. There exist different metrics whicli 
measure the centrality of a node (among others: degree, betweenness, closeness, eigenvector) fS], each 
one highlighting different facets of this generic concept. 

Among the empirical analyses, there is a body of literature showing that centrality explains the 
role, importance, or payoffs of the agents constituting the network: in informal structures within 
organizations, the importance of people is related to their betweenness centrality [5]; the outcome 
in education tests is related to the betweenness and closeness centralities of the individuals in the 
network [SQ From the theoretical side, in Ref. [T it was shown that, in a game where individuals have 
to maximize their utility - and under mild assumptions of concavity of the payoff function -, their 
gains result proportional to their Bonacich centrality |4] . Because of this, if individuals can alter their 
neighborhood, the myopic best response strategy is simply to connect to the neighbor who increases 
their centrality the most. 

Interestingly, in Ref. [TT] it was shown that when individuals strive to be as central as possible, the 
exact measure of centrality is irrelevant, and the dynamics yields a network which has the property of 
nestedness: the neighborhood of any node contains the neighborhood of the nodes which have a lower 
degree. In this kind of networks, the ranking of nodes according to their centrality is the same, regardless 
of the centrality measure considered 10 . Because of this property, the network dynamics is independent 
of the exact centrality measure considered. Interestingly, nested structures have been found in inter- 
organizational networks of research and development (R&D) alliances [21lll7j . in interbank payment 
networks [19^ and in firm competition under oligopolies [S]. This kind of structures will be precisely 
the ones emerging in the social climbing game. In this respect, our results confer stability to those of 
Ref. [To] and generalize them in non-trivial ways. 

This study allows us to point out that the assumption that individual freedom promotes social 
mobility is a non-trivial one. This is because the structure of a society, while constraining the set of 
opportunities that are available to individuals, depends on the very incentives of individuals in complex 
ways. In this paper we show that this interplay may produce very "rigid" societies, with extremely 
low social mobility, characterized by persistent inequalities between a priori equal individuals. The 
understanding of this phenomenon hinges on the concept of ergodicity breaking that occurs in strongly 
interacting systems, when a symmetry - here related to the a priori equality among individuals - is 
spontaneously broken. This phenomenon, well known in statistical physics, is an emergent collective 
property, and it manifests only when the system is large enough. Remarkably, we find that persistent 
inequality with low mobility occurs precisely when the quest for "power" - i.e. for occupying the most 
central or important place in the social hierarchy - becomes a dominant component of what motivates 
the behavior of individuals. 

In words, our model epitomizes an apparent positive feedback between the intensity of the efforts 
of individuals to "climb" the social hierarchy and the structure of society: on the one hand, the more 
a society is hierarchically structured, the easier it is for individuals to understand how to climb it. On 
the other, the efforts of agents to climb the hierarchy reinforce the social ranking as individuals rewire 
their links from less to more influential individuals. We discuss this interplay in a highly stylized model 
of a society, that while being very far from realistic, serves as a proof of concept and allows us to unveil 
the mechanism responsible for the emergence of a persistent inequality in a transparent manner. 

In addition, this approach makes it possible to deploy techniques used in statistical mechanics [ini 
[T5] in the context of network properties. Similar models to the one considered here have been discussed 
in Refs. [1311201) that, however, focus mostly on topological properties of the emerging networks. 

The rest of the paper is organized as follows: In Sect. ^ we introduce the model and discuss the 
main properties of the dynamics and its associated a global potential function; related to these results, 
in Appendix [Al the ergodicity of the system is proved. Later, in Sect. [3] we show the results of extensive 
numerical simulations that portray the characteristic behavior of the system. Finally, in Sect. [4j the 
conclusions are drawn. 



2 The model 

We consider a system composed of N individuals, who are connected through a network which consists 
of exactly M links. The network is undirected and thus can be specified in terms of a symmetric 



^ A more comprehensive list can be found in Ref. 
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adjacency matrix a = with elements Qij = Qji = 1, if i and j are connected, a.y = aji = 

otherwise. Agents receive opportunities to use their hnks in order to get in contact with more 
"influential" members of the society, in brief to climb the social network. 

As a measure of importance of the individuals, we take the number of his/her partners, [^i.e. the 
degree ki = ^jCiij- Therefore, a measure of the "social capital" of agent i is given by the following 
local utility function 

N N N 

Ui = ^ aijttji + ^ an = ^ aijk-j + ^h, (1) 
e=i j=i 

whose first term counts the total number of friends of his/her neighbors, while the second one is simply 
his/her degree and /i is the relative weight between these two term^ The efforts of agents to climb 
the social hierarchy can then be formalized in the maximization of the utility Ui. 
We then define the dynamics as follows, 

1. At any time, an agent i is picked at random together with one of her neighbors, j. Then, a neighbor 
of j is selected at random, £ ^ i. 

2. If £ is already connected to i, nothing happens. Otherwise, with probability 

P{i,3)^{^,l) = ^ ^ , (2) 

where Aui is the change in i's utility, the link is replaced with another one (i,^). 

The first step models random encounters between agents through their network of interactions. In such 
an encounter, agent i gets to know a friend £ of j, as well as his/her importance (the number of Ps 
friends). At this point, if the only concern of an agent were his/her social status, he/she would always 
choose the move implying an increase in his/her local utility. However, in reality, the social network 
supports many more functions and benefits than those that are directly related with social status. This 
can be accounted for by a random utility model where only part of the utility is observecj^ 

In the second step, the randomness of the utility is encoded in the so-called intensity of choice 
parameter _/?, which reflects the prevalence of the quest for social status in how agents decide about 
their linka^ In particular, in the limit f3 — >■ oo, a move implying a decrease in the utility function is 
never accepted. This means that the social status is valued so highly by the agents that everything else 
is unimportant. On the contrary, for /3 = the probability of accepting a move implying a decrease of 
the utility function is 1/2, meaning that the social status has negligible importance with respect to the 
unobserved part of the utility. The general question addressed is then how strong should the intensity 
of choice f3 be in order for a social hierarchy to form and be maintained in the long run? 

It is worth to remark that if the utility of agent i increases when rewiring the link to (i,^), 
then the utility of agent j decreases, while that of agent £ increases. This embodies the fact that 
the formation of a new link needs the consent of both parties, but their removal can be unilateral. 
Moreover, this also implies that the utility of an agent is also increased by the decisions of others. 
In their quest to become central, the agents increase the likelihood to be selected by others as new 
partners. 

Notice finally that the number of links is conserved in the dynamics. Hence the density of links is 
the second important dimension that we shall explore, in order to understand how the structure of 
social organization depends on it. 

^ Other measures of centrality can be taken but, as observed in Ref. [10], these rank individuals in the 
same order in strongly hierarchical networks, that will be stable over time as we shall see later. Conversely, 
unstructured networks correspond to random rankings with no stable order, with respect to all centrality 
measures. 

* As will be clear in the following, the second term in (|T]) is irrelevant for the dynamics, but not for the 
interpretation of the local utility. For example, consider the limit case in which one has a central node connected 
to — 1 nodes whose only neighbor is precisely the central node. In this case the first term in (IT]) is equal to 
N — 1 both for the central node and for the remaining ones. One way to remove this degeneracy islo introduce 
a term proportional to the degree. 

^ See [7j p. 33 for an explicit derivation of the stochastic choice model from random utility models. 

® For example, Adam may be reluctant to interrupt his relation with Bob, despite his low rank in society, 
because he is his only friend who shares his interest in Japanese paintings. 
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2.1 Properties of the dynamics and potential function 

There are some remarkable features of the dynamics of the model introduced in the previous sections. 
We detail them now. 

First, it is easy to see that the dynamics introduced preserves connected components. Indeed, 
nodes are never disconnected by the dynamics because, even if they have just one link, this will not 
be rewired because the neighbor upstream has no second neighbor where to rewire. Therefore, without 
loss of generality, we restrict attention to the case where M > N — 1 and the network is composed of 
a single connected component. Networks composed of disjoint components remains disjoint under the 
dynamics above, hence their dynamics can be considered independently. 

As a consequence of this fact, and since both the number of N nodes and M edges is left unchanged 
during the evolution of the system, the number of fundamental cycles in the graph is conserved. This 
is a simple corollary of the proof in [31 Ch. 2], where it is shown that the number of fundamental 
cycles in a graph is equal to M — N + K, where K is the number of connected components. 

Interestingly, the dynamics of the model admits a potential which is just the global utility, i.e. the 
sum of the utilities U = ^ ^ Ui . Indeed, let us consider the change Au^ in the utility of the agent x 
when the rewiring into (i,£) occurs. Depending on the position of x in the network, the following 
changes are obtained: 



Aui ~ ki — kj + I (3a) 

Auj ^ 1 — ki — fi (3b) 

Aue = fci — 1 + /i (3c) 

Auh = -l yhedj\{i,e} (3d) 

Aug = +1 yged£\ {]} (3e) 

Au^ = +0 yx^i,j,e,x (^djUdi, (3f) 



where dx is the set of the neighbors of x, before the move. 

In the total variation of the utility AU = Aux, the term Auh appears kj — 2 times, while the 
term Aug appears ki — 1 times, because k^ is the degree of the node x before the rewiring. Gathering 
all the contributions one has: 

AU = Au^ + Auj + Aug + {kj - 2)Auh + - l)Aug 

= 2{ki - kj + 1) = 2Aui . ^ ' 

The last point implies that the dynamics should converge to the equilibria of a system with Hamiltonian 

i i 

and fixed density of links at temperature 2//Q Notice that the second term does not play any role, 
being ki a fixed quantity in our case. Indeed the dynamics in Eq. ^ samples the Gibbs distribution 
P{d} cx e"^*^^"-*/^, provided the dynamics is ergodic, which is proven in Appendix [a| As a consequence, 
our model is closely related to the 2-star model [T6 1 IT5 ] , that in the case where the density of links is 
not fixed, exhibits a sharp phase transition. This result suggests that there might be a phase transition 
also in the model we study in this paper. It also provides a microeconomic derivation for the 2-star 
model, that is widely studied in sociology. 

3 Numerical simulations 

In order to investigate the behavior of the model, we performed extensive numerical simulations sam- 
pling the Gibbs distribution P{a} oc e"^'^^"-'/^ using the Metropolis algorithm based on the rewiring 
moves introduced in Sect. [2] All the results to be presented throughout the rest of this section were 
obtained, for each value of /3, by performing R rewiring proposals per node, and we checked that the 



The factor 2 comes from the fact that the variation of the global utility is the double of the variation of 
the local utility. 



Fig. 1 Snapshot of networks of the social climbing game for — 100, A/ — 125 for /3 = 0.03 (left panel) and 
/3 = 0.1 (right panel). Size of the nodes is proportional to the degree. 



value R = 5 - 10^ is large enough to always ensure the attainment of an equilibrium state. Fig. [T] shows 
two typical realizations of the social network for small and large values of (3 (see caption for more 
details). Fig. [I] suggests that, as anticipated in the previous section, the social climbing model under- 
goes transition from hierarchical to random structures. In the following, we will show the presence of 
a phase transition between these two states. 

In Fig. [2] we show the largest degree of the network ^ — max^ ki as a function of the inverse 
temperature /? for systems with N = 100 nodes and M = 110, 200, 300, 500 links. As can be seen, in 
all cases the system actually undergoes a transition, going from a phase where the largest degree ^ is 
roughly of order 1 — 10 (depending on the relative size of N and M) to a phase where the largest degree 
is of order N. These observations qualitatively match the findings of [13], where a prediction for the 
critical temperature Tc — 1/ Pc characterizing this phase transition was also derived, from combinatorial 
arguments, for networks with average degree k = 2M /N < 2, i.e. for disconnected graphs. The nature 
of the phase transition depicted in Fig. [2] is further investigated in Fig. [3j where the relation between 
the average and maximum degree (divided by N) is explored for a network with N = 100 nodes at 
a fixed inverse temperature /3 = 0.03. As can be seen, a coexistence (in statistical sense) between the 
disordered and ordered phases is observed, signaling a first-order phase transition. 

Following |T3j, let us define a star as a node whose degree is of order N. Then, it is immediate to 
figure out that, depending on the ratio M/N, different number of stars might emerge in the network 
for temperatures lower than Tc. Clearly, in the case N = 100, M = 110 (i.e. k only slightly larger 
than 2), the appearance of a star {<P ~ 100 in this case) below the critical temperature leaves very 
few links to be distributed amongst the remaining nodes. On the other hand, increasing the number 
of links provides enough room for the emergence of a larger number of stars. In other words, it is 
intuitively reasonable to expect a system with an average degree fc ~ 2n to produce, for sufficiently low 
temperatures, exactly n stars. In order to support such an intuitive line of reasoning, we computed the 



inverse participation ratios (IPRs) of the degree sequences v = (fci, fc2, • ■ • , ^a')/v of several 
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Fig. 2 Dependence of the largest degree <l> (divided by TV) in the social climbing network as a function of the 
inverse temperature (or intensity of choice parameter) /3. The different curves refer to = 100 and M = 110, 
200, 300, 500. For each value of /? the reported values of <l> are obtained by averaging over 100 networks. An 
abrupt change in is observed in all curves after a threshold value of P, with <1>/N going from low values to 
values close to one, signaling the emergence of a star, i.e. a link with 0{N) links, in the network. 
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k 

Fig. 3 Density plot depicting the relation between the average and maximum degree (divided by A*') at inverse 
temperature P — 0.03 and number of nodes A'^ = 100. High (low) values are in red (blue). Clearly, for low values 
of k the network is in the disordered phase, while for high values of k it is in the ordered phase. Coexistence is 
observed between the ordered phase, characterized by a largest degree of order N, and the disordered phase. 



networks with different numbers of nodes and links. Given a normalized vector v, its IPR is defined as 



The IPR of a completely localized vector, say v — (1,0, ... ,0), is equal to one. On the other hand, 
the IPR of a fully delocalized vector, whose components are all equal to Vi ~ is of order N. 
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Fig. 4 Inverse participation ratio of tlie normalized degree sequence /(v) as a function of the inverse tempera- 
ture /3. Different curves refer to networks with A'^ — 100 and M — 110, 200, 300. Results obtained by averaging 
over 100 networks. For large enough /3 one recovers the maximal number of stars. 



In our case /(v) gives an estimate of the number of dominant nodes in the network. In Fig. |4]we plot 
the IPRs /(v), as functions of /3, for N = 100 nodes and M = 110, 200, 300. As can be seen, the 
IPR of the sparsest network, i.e. the one with M = 110, essentially drops down to one right below its 
critical temperature. On the other hand, systems with a larger number of links undergo a less trivial 
evolution: after the initial drop below the critical temperature, the IPR increases and eventually reaches 
a steady state. In the example shown in Fig. |4] the system with M = 200 links reaches a steady value 
/(v) ~ 2.12 ± 0.03, whereas the system with M = 300 reaches /(v) ~ 3.28 ± 0.06 (where the errors 
represent the 68% confidence intervals obtained by averaging over 100 networks), and such values 
clearly show that the maximal number of stars allowed by the relative sizes of TV and M has been 
achieved. Moreover, these observations are consistent with the small temporary decrease of the largest 
degree <P which can be observed in Fig. [2] for systems with k > 2 when the inverse temperature is 
slightly larger than its critical value. 



3.1 Correlations and social mobility 

As already explained in Sect. [TJ one of the goals of the present paper is to model the positive feedback 
mechanism between the individuals' effort to climb the social hierarchy and the subsequent reinforce- 
ment of the social hierarchy itself. Suppose that a given social network reaches its equilibrium state, 
at a certain inverse temperature after steps of the social climbing dynamics described in Sect.[2j 
Let us denote the corresponding graph's adjacency matrix as a{to). Then, one way of quantitatively 
describing how mobile or "frozen" a society is would be to assess the level of correlation, according 
to some proper notion, between a{to) and a following configurations a{t), where t = to + At for some 
positive At. We will now measure correlations by means of Kendall's rank correlation coefficient. Given 
the joint set of all matrix entries in a(to) and a(t), let us focus, for example, on entries and {h,i) 
in both matrices. Then, if both aij{to) > ahe{to) and aij(t) > ahe{t), or if both aijit^) < ahi{to) 
and aij{t) < Uhiit), the pairs (a^j (<o), aM(io)) and {aij{t),ahi(t)) are said to be concordant. On the 
contrary, if ajj (to) ^ a/i£(io) and a„ (t) ^ aM{t) the pairs {aij{to),ahe.{to)) and {a,j{t), aneit)) are said 
to be discordant. Of course, since the adjacency matrix entries equal zero or one at each time, ties will 
often happen either at time to or at time t (or at both times). Kendall's correlation coefficient r reads 

^i^t) - (6) 
^C + D + Tt„^C + D + Tt' 
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Fig. 5 Kendall's r coefficient (see ml) measurements for networks with TV = 100 and M = 300. All measure- 
ments are performed between an inmal equilibrium configuration a{to) and later configurations a{to + nAt), 
with to — At = 5 ■ 10® Monte Carlo steps. The different curves refer to inverse temperatures /3 = 0.02, 0.04, 
0.06, respectively corresponding to values below, slightly above and well above the critical value for the system 
under study (see also Fig. [2]|. Results obtained by averaging over 100 networks. 



where C {D) is the numbers of concordant (discordant) pairs, whereas Ttg (Tt) denotes the number of 
time-to (time-i) ties. Pairs where ties happen both at to and t are not taken into account. 

In Fig. [5] a few examples of Kendall's r coefficient's time evolution are sketched. All plots refer 
to networks with N = 100 nodes and AI — 300 links. Here, At = to = 5 ■ 10^ elementary Monte 
Carlo moves, i.e. rewiring proposals. As can be seen, when the social climbing game takes place for 
temperatures higher than the critical one, Kendall's r quickly starts to fluctuate around zero, denoting 
no genuine correlation between configurations distant (in time), denoting a large social mobility. On 
the other hand, for temperatures slightly lower than the critical one, Kendall's t remains significantly 
larger than zero over several time lags. However, a downward trend is clearly visible in this case, 
meaning that for temperatures T < Tc social mobility is recovered after a sufficiently long time. On 
the contrary, for temperatures significantly lower than the critical one Kendall's r essentially remains 
constant and very large (i.e. close to one) over large time lags, hinting at an extremely reduced social 
mobility, possibly preventing the majority of individuals from climbing the social ladder. 

The above considerations on individuals' mobility in the social climbing game can be further clarified 
and understood more deeply. For these purposes, let us denote as qi the fraction of agents who, at a 
given time, have a strictly lower degree than agent i, i.e. 

1^-^T.(^ik^~k,), (7) 

where 9 denotes Heaviside's step function (whose value in zero is here fixed as 0{O) — 0). The variable 
defined in ([7| clearly represents a suitable definition of the social ranking, hence the social status, of a 
given individual in the network. Thus, a reasonable measure of the individuals' mobility in the social 
climbing game is given by the change in the quantity defined above over a certain time lag At, i.e. 
Aq,{At) = q^{t + At) - q,{t), for i = 1, . . . , TV. 

In Fig. [I some typical behaviors of the q index defined in ([7| are shown. All examples refer to 
networks with N = 100 nodes and M — 1000 links. In such plots, the average of the change Aq is 
shown as a function of g, in order to provide information about the typical social mobility over a time 
lag At for an agent whose social ranking at the beginning of such a time lag is quantified by q. As can 
be seen, depending on the preference for social status, i.e. on the inverse temperature parameter /3, very 
different situations can happen. In a rather disordered society (low values of /3) the relation between q 
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Fig. 6 Relation between the q index defined in n\ and its variation Aq over a given time lag At = IC* Monte 
Carlo steps for a network with A*' — 100 and Kl — 1000. The different curves refer to inverse temperatures 
/3 = 0.02, 0.08, 0.16, respectively corresponding to values below, slightly above and well above the critical 
value for the system under study. Points refer to the average variation Aq over an equally spaced grid of q 
values (going from to 1 in steps of 0.05). Results obtained by averaging over 100 networks. Shaded area (for 
/3 = 0.02) and error bars (for P = 0.08, 0.16) represent the central 68% of the events. Points and error bars 
relative to different values of j3 have been shifted to enhance readability. 



and Aq is clearly linear, and does not depend strongly on the time lag size At. In particular, it can be 
seen that, on average, individuals sitting at the bottom of the ranking typically end up higher in the 
social ladder after some time, whereas individuals sitting atop the hierarchy are prevented from keeping 
their social status intact for a long time. When the preference for social status crosses its critical value, 
such a picture starts changing quite dramatically. For values of /3 slightly larger than the critical value 
Pc agents with low degrees still have a chance to climb up the social ladder, especially over rather long 
time lags, whereas the dominant individuals (g > 0.9) typically get to keep their social ranking. Also, 
as can be seen from the right plot in Fig. |6j when the critical threshold is crossed the social network 
becomes "fragmented", as the q index is no longer defined over the whole interval [0, 1]. In a strongly 
ordered society, i.e. (3 well above its critical value, agents with low degrees are almost completely stuck, 
and all of the social mobility happens in the top half of the social network, i.e. amongst agents with 
q > 0.5. These results complement, at a "microscopic" level, those presented in Fig. [5j 



4 Conclusions 

In summary, we have discussed a very simple model for the dynamics of a social network where the 
agents' quest for high status in the social hierarchy reinforces the latter while reducing social mobility. 
The model is very stylized and far from a realistic description of social dynamics. Yet, it captures 
some key ingredients that are enough to reproduce stylized facts known at least since the work of 
Vilfredo Pareto [H] . Namely, Pareto observed that societies tend to organize in a hierarchical manner, 
with the emergence of "social elites"]^ Our model, as well as Refs. jTTlITU] , provides a formal framework 
showing that individual incentives for high social status are enough to confer this property to the social 
network, even in the absence of explicit discrimination of particular groups (e.g. cast system or racial 
segregation) or preferential biases (e.g. hereditary rules). In addition, we find that the hierarchical state 
is remarkably stable, with suppressed social mobility in the upper and lower parts of the hierarchy. 
Notably, Pareto himself observed that social mobility is higher in the middle classes [Uj. Our model 

® A similar concept of "power ehtes" has been discussed in [T5] . 
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neglects important dimensions, such as wealth or political power that, however, likely contribute to 
reinforce our results. 

Secondly, we show that the social climbing game admits a potential function, thereby allowing us 
to deploy techniques and concepts of statistical mechanics to understand the behavior of the system. 
Statistical mechanics provides a natural language for discussing collective properties of societies. For 
example, the emergence of a social hierarchy in a system of a priori identical individuals is an example 
of spontaneous symmetry breaking, whereby the associated loss of ergodicity accounts for the reduced 
social mobility. 

The present paper was mostly focused on investigating the model via numerical simulations. The 
mean field approach discussed in Ref. [TOlfTS] is not applicable in our case, because the density of links 
that plays the role of an order parameter in Ref. P^[T5] is fixed in our case. Indeed, the phenomenology 
we find is different from that of Ref. [16.. 15] as we do not find evidence of hysteresis phenomena: there 
is no range of parameters where the disordered and the ordered societies are both stable. We speculate 
that this might be related to the fact that in the social climbing game there are mechanisms by which 
a social hierarchy can "nucleate" gradually in an ordered society, by forming social elites that grow 
over time. 

This and other issues can in principle be addressed within more sophisticated statistical mechanics 
approaches. In this respect, it is worth to mention that it is possible to map the problem into that 
of an interacting lattice gas that possibly admits for a full and exact statistical mechanics treatment. 
Work in this direction is currently in progress. 
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A Ergodicity of the dynamics 

Let r^(iV,M) be the space of connected graphs with A'^ vertices and M edges. In order to prove ergodicity, 
we have to show that, with a finite number of basic moves, we can reach any connected graph in (N, M), 
starting from another arbitrary graph in the same set. Before delving into the technical details, we give a simple 
intuitive sketch of this proof. 

For a finite value of /3, the dynamics consists of reversible moves as the one depicted in fig. (TI such moves 
can be thought of as a "sliding" of the edge Cik on the path of length one from vertex k to vertex j. 

The key observation to prove the ergodicity by induction is that, since the graph is finite and connected, there 
always exists a path of minimum length that connects two arbitrarily chosen vertices in the graph. Then, we 
can proceed in three steps. 

1. Let there be two graphs in F'^ {N, M) which differ from each other only by an edge incident on the same 
vertex, Vk- We first prove that by means of basic moves, we can transform one into the other. To do so, it 
sufiices to slide the edge along the path that connects the other end of the edge, which we know to exist 
because the graph is connected (Prop. 1, Prop. 2). 

2. Let there be two graphs in F'^{N, M) that differ by an edge with arbitrary ends. By applying the previous 
step twice, we show that there exists a finite set of moves that allows us to reach one configuration starting 
from the other (Prop. 3). 

3. Finally, let there be two arbitrary graphs in F'~' {N, M). Moving one edge at time, we show by induction 
that it is possible to reach one graph starting from the other with a finite number of moves. Thus, the 
ergodicity is proved (Prop. 4). 

We now proceed with the detailed Proof. 

Definition 1 (c-swap) Let us choose a labeling for the space of vertices V = {vi, . . . ,vn} and an induced 
labeling for the edges E = {eij} where eij = eji = {vi,Vj) denotes the undirected edge between Vi and Vj. Let 
us define a transformation a-j : F'^'{N,M) i— >■ F^'{N,M), called corner swaps (c-swaps), as following 



g'^iv,E') 



(8) 



such that 



E' = 





otherwise 



(9) 
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Proposition 1 Let Q = {V, E) and Q' = {V, E') be two graphs in r'^{N,M) that differ by an edge incident 
on the same vertex, i.e. \E\ = \E'\ = M, \E' n E\ = M ~ 1, E\ E' ^ {eik} and E' \ E = {cij}, and such that 
the shortest path P from vt to Vj does not contain neither Vi nor any of its neighbors. 
There exists an integer I and a finite sequence of graphs in F'^ {N, M), Q" such that: 

(i) g = and g' = gK 

(li) For all < n < I there exist adjacent vertices Ufc„,Ufe„+i such that g"'^^ — cr*^"^^ (C/"), where ko — k and 
ki =j. 

Proof Let / be the length of P. 

Let Vki be the unique neig hbor of Vk that lies in P. If we set g^ = al^^ (g). The c-swap reduces the distance 
between Vi and Vj , since the neighbor of Vk that lies in P must have a distance I — 1 from Vj . We reiterate the 
procedure on g^ and obtain in such a way a sequence of graphs that satisfies property (ii). Now, since at any 
step the length of P diminishes by 1, after the l-th step, in the graph g' Vi and vj will be neighbors. Thus, 
since no other edge was changed by applying c-swaps, g' = g' proving property (i). 

Proposition 2 Let g — {V, E) and g' = {V, E') be two graphs in r^'{N,M) which differ by an edge incident 
on the same vertex, i.e. \E\ = \E'\ = M, \E' n E\ ^ M - 1, E \ E' ^ {e,fc} and E' \ E ^ {e,j}. 
There exists an integer I and a finite sequence of graphs in F'~^ (N, M), g" such that: 

(i) g ^g° and g' = 

(li) For all < n < I there exist adjacent vertices Vk„,Vk„+i such that C/""*"^ — o-l^"^^{g"). 

Proof Let P be the shortest path in g from Vk to Vj that does not contain (vk, Vi). 
There are four possible cases : 

(i) P does not contain neither Vi nor any of its neighbors other than v^. The thesis is proven applying 
proposition [T] directly to P. 

(ii) P contains Vi. Let Pi be the shortest path from Vk to Vi that does not contain the edge {vk,Vi). let 
P2 be the shortest path from Vi to Vj, clearly P = Pi © P2, where © is the path concatenation. Since 
by construction there are no neighbors of Vk in P2 (otherwise P would not contain Vi) we can apply 
proposition [1] and reach g" = {V, {E \ {cki}) U {ekj})', on the other hand there cannot be neighbors of v_j_ 
in Pi (otherwise there would be a shortest path not containing Vi) and thus applying again proposition [l] 
along Pi we reach g' proving the thesis. 

(iii) P does not contain Vi but two of its neighbors, c and / such that c ^ Vk, f ^ Vk and |c, Wfe| < |/, Wfe|, where 
]-, - I represents the graph distance between two vertices. We first note that c and / must be neighbors, 
otherwise P should include Vi. Then, as in case (ii), by minimality we can write P = Pi © (c, /) © P2 
where Pi is the shortest path from Vk to c and P2 is the shortest path from / to Vj. It is easy to see that 
Q2 = {vi, /)©P2 is a shortest path from Vi to Vj-. if it were not so, there would exist a path Q'2 from Vi to Vj 
strictly shorter than Q2, but in that case Pi©(c, Vi)(BQ'2 would be a shortest path from Vk to Vj containing 
Vi, in contradiction with our hypotheses. A similar argument holds for Q\. As before, since, by minimality, 
there cannot be neighbors of Vk in P2 it is possible to reach the graph g" — {V, {E \ {eki}) U {ekj}) by 
applying proposition [T] to Q2; since by minimality there cannot be neighbors of Vj in Pi, we can apply 
proposition [1] to g" along Q2 and reach g' proving the thesis. 

(iv) The shortest path P contains only one neighbor of Vi other than Vk, let us call it m. As before, P — Pi ©P2 
where Pi is the shortest path from Vk to m and P2 is the shortest path from m to Vj . Since by construction 
there cannot be other neighbors of i in P2 , we can apply proposition [l] to P2 and reach the graph 
g* — (y, (£ \ {eim}) U {cij}. On the other hand, by construction there cannot be neighbors of Vi in Pi 
other than Vk and thus we can apply proposition [l] to P2 and reach g' proving the thesis. 

Proposition 3 Let g = {V, E) and g' = {V,E') be two graphs in F^{N,M) such that \E\ = \E'\ = M and 
|P n P'l = M - 1. Let us assume that, in particular, E = {e^j} U (P n E') and E' = {e^fc} U {E n E'). 
Thus there exists an integer I and a finite sequence of graphs in P^ {N, M), g" such that: 
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(i) g = and g' = g^ 

(ii) For all < n < I there exist adjacent vertices Vk„,Vk„+i such that g"'^^ = cr^^"^^ (C/"). 

Proof Let us define the graph g" — {V,E") such that E" = {E \ {cij}) U {cih}- Applying proposition [2] first to 
graphs g and g" and then to graph g" and g' proves the thesis. 

Definition 2 (5-swap) Let g = {V, E) and g' = {V,E') be two graphs in r^(N,M) which differ at most 
by an edge, that is such that \E\ = \E'\ = M and \E n E'\ = M — 1. Let us assume that, in particular, 
E = {cij} U{En E') and E' = {enk} U n E'). 

We define a global swap or g-swap of the edge Cij to the edge ehk a transformation such that: 

e;' = rf/(e;) (lo) 

Proposition [3] simply states that any global swap can be obtained as the composition of a minimal set of corner 
swaps between adjacent vertices. 

Proposition 4 Let g — {V,E) and g = (V, E') be two graphs in r^(N,M). There exists an integer d and a 
sequence of graphs g"{V,En) in r'~^'{N,M) such that: 

(i) g = g° and g' = g'' 

(ii) For all < n < d there exist four vertices Vi , Vj ,Vh and Vk such that C/""*"^ — I]^^{g'^) 

Proof Let Z = (V^, Z — E n E'), and let us define S = \Z\. We proceed by induction on the number S. 
Base case If 5 = A/ — 1, the Thesis is trivially true because of Proposition [3] 

Inductive step Let us assume that the Thesis holds for 5 = M — d, we want to show that this implies that it 
also holds for 5 = M — d—1, with d < M — 1. Let us assume that g = {V, E) and g' = {V,E') are such that 
|£"n£;| = M-d-l. Let e^- G E\{Er]E') and Chk G E'\{Er)E'). Moreover, let E" = {E\{e^j})U {ehk}- 
By construction, \E n E"\ = M - 1 and \E' n E"\ = M - d. Finally, let g" = {V,E"). Since g" and 
g' differ by M — d edges, by inductive assumption there exists a sequence C/*, with i £ [0,d], such that 
5" = g' and g'^ = g" , that satisfies the Thesis. Moreover, by Proposition [3] there exists a p-swap such that 
g = E'^{{g"). Thus, the complete sequence 5' = 5°, e\ ■ ■ • , G" = Q'^,Q = G''^'^ satisfies the Thesis. 

Propositionjl] and proposition [3] state simply that any two connected graphs with the same number of edges 
can be obtamed one from the other applying a finite sequence of c-swaps. Moreover, since the number of edges 
is finite, then there must be a minimal sequence of c-swaps that connects any two of such graphs. Since for 
finite /3 model all c-swaps are allowed with certain probability this proves the ergodicity. □ 
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